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Abstract 

Virtual radiative corrections due to the long range Coulomb forces of heavy nuclei with charge Z 
may lead to sizeable corrections to the Born cross section usually used for lepton-nucleus scattering 
processes. An introduction and presentation of the most important issues of the eikonal approxima- 
tion is given. We present calculations for forward electroproduction of rho mesons in a framework 
suggested by the VDM (vector dominance model), using the eikonal approximation. It turns out 
that Coulomb corrections may become relatively large. Some minor errors in the literature are 
corrected. 
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1 Introduction 



Due to the small size of the elemental electric charge e = -v/ivra, it is for most electromagnetic 
elementary particle reactions fully sufficient to calculate only in Born approximation. But this is no 
longer true when heavy nuclei are involved, like e.g. lead, where the relevant perturbation expansion 
parameter aZ ~ 0.6 is not small. It is then possible that even at high energies, the ratio of the exact 
and the Born cross section does not approach unity. One prominent example for this observation is 
the photoelectric effect, where a single photon knocks out an electron from the K shell, or electron 
positron pair production by a single photon incident on a nucleus, where the electron is captured into 
the K shell |2j. For small aZ and high photon energy e^, the cross section for these processes is 
given by the Sauter formula [S] 

do ~ Z^A^ — , (1) 

where Ac is the electron Compton wavelength. Applying the Born approximation in the usual sense, 
i.e. by the use of plane waves for the positron, was attempted by Hall and Oppenheimer pP| already in 
1931 in the case of the photoelectric effect. But a central difficulty in the treatment of the photoeffect 
arose from the distorted wave function of the ejected electron. The bound state wave function depends 
on aZ and only terms of relative order aZ survive in the Born matrix element. Terms of this order 
also come from the continuum wave function. Using exact wave functions for the continuum electron, 
it turns out that the cross section calculated from plane waves for the ejected electron is already wrong 
for small aZ. It is therefore not astonishing that Coulomb corrections can become very large even at 
high energies. For the derivation of the Sauter formula, nonrelativistic bound state wave functions and 
approximate Sommerfeld-Maue wave functions for the ejected electron were used. Boyer |5j improved 
the calculation by using exact bound state wave functions and Sommerfeld-Maue wave functions for 
the continuum wave function. The calculational advantage of using Sommerfeld-Maue wave functions 
relies in the fact that the exact continuum wave functions are only available as an expansion in partial 
waves, and for higher energies summation over a large number of terms is necessary in order to compute 
exact matrix elements. But Sommerfeld-Maue wave functions are valid for the Coulomb potential of 
point-like charges, and hence do not take into account the finite size of the nucleus. We therefore 
pursue a similar strategy, namely the eikonal distorted wave approximation, where the advantage is 
that also finite size effects can be taken into account. This is clearly necessary for nuclear processes, 
whereas the relevant physical length scale of the photoelectric effect is given by ~ {aZ)~^Xc and much 
larger than a typical nuclear radius. 

In the following part of this paper, we revisit the eikonal approximation by giving a condensed 
introduction to the subject, where we also point out the most important facts related to possible 
improvements concerning the range of validity of the method in its basic form. In sect. 3, we apply 
our calculation to coherent vector meson electroproduction. It turns out that the Coulomb corrections 
are indeed large for heavy nuclei, and that it is necessary to use a rather accurate description of the 
electrostatic nuclear potential. We correct some errors in the literature and rederive the eikonal result 
given in 0. 

2 Eikonal Approximation 

For highly relativistic particles in a potential V with asymptotic momentum p, one may neglect the 
mass of the particle y|, \p \ ^ m), such that the energy-momentum relation reduces to (h = c = 1) 

I^E -Vf =p'^ + rr? ^ E -V = p,p=\p\. (2) 
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The classical relation for the (initial) momentum of the particle pi = piPi can be taken into 
account approximately in quantum theory by modifying the plane wave describing the initial state of 
the particle by the eikonal phase xi(0 



Xiir) = - / V{f+pis)ds. (4) 



-oo 



If we choose the initial momentum parallel to the z-axis pi = p^e^, the phase is 

z 

Xi(r1 = - / V{x,y,z')dz'. (5) 



The z-component of the momentum then becomes 

- id,e'P-'+'^' = (pI - V)e'P>+'^' , (6) 

and further analysis shows that also the transverse momentum modification is well approximated. 
Analogously, for the final state wave function follows 

where 

oo 

X2{^ = - jv{r+Pfs')ds' . (8) 


For the sake of simplicity, we consider spinless particles first. The spatial part of the free charged 
particle current 

- ie[e-'Ps''Ve'P^'" - (Ve-^^^^e^P^n = <Pi + Pf)e'^^'-^i^'' (9) 

becomes 

e{p, +pf + Vxi - Vx2)e*(P^^-^^^')^^+''^('^ , (10) 

where e is the charge of the particle and x(^) = Xi(^ + X2(^)- The spatial part of the current now 
contains the additional eikonal phase, and the prefactor contains gradient terms of the eikonal phase 
which represent essentially the change of the electron momentum due to the attraction of the electron 
by the nucleus. 

So far we have only considered the modification of the phase of the wave function, and for many 
problems, like the one treated in this paper, this is a sufficient approximation. It has been applied to 
elastic high energy scattering of Dirac particles in an early paper by Baker jSj. However, the method 
leads, e.g. for quasielastic nucleon knockout scattering of electrons on lead with initial energy ei = 300 
MeV and energy transfer uj = 100 MeV, to errors up to 50% in the calculated cross sections. The 
reason is that also the amplitude of the incoming and outgoing particle wave functions is changed 
due to the Coulomb attraction. This fact is related to the classical observation that an ensemble of 
negatively charged test particles approaching a nucleus is focused due to its attractive potential. For 
the sake of completeness, we give here a short discussion of this effect. 
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Knoll [2 derived the focusing from a high energy partial wave expansion, following previous results 
given by Lenz and Rosenfelder jlUl lllj . E.g., for the incoming particle wave in the vicinity of the 
nucleus he obtained an approximation for the electron wave function (up is the constant free electron 
spinor) 



ipi = e'^+{p' /p)e'P''^{l + air^ - 2a2p'r + iair^p'f+ ia2[{p' x f)^ + a{p' x r)]}up, (11) 

where a describes spin changing effects. A similar formula holds for the final state wave function. The 
so-called effective momentum p' is parallel to p and is given by the classical momentum of the electron 
in the center of the nucleus, i.e. p' = \E — V{0)\. The parameters ai^2 depend on the shape of the 
nuclear electrostatic potential. For a homogeneously charged sphere with radius Rg they are given 
by ai = —aZ/6p'Rg, 02 = — 3aZ/4p'^i?^. The increase of the amplitude of the wave while passing 
through the nucleus is given by the —2a2p'r-teicm. The air^-Term accounts for the decrease of the 
focusing also in transverse direction. Imaginary terms like ia2{p' x r)^ describe the deformation of the 
wave front near the center of the nucleus. They could be obtained correspondingly by an expansion 
of the eikonal phase in that region, and 5+ is the eikonal phase in the center of the nucleus x(0)- 

Apart from the spin structure which is absent for scalar particles, the eikonal approximation relies 
basically on the same strategy for particles with or without spin, namely on the modification of the 
wave function by the eikonal phase. But there are also gradient terms present in the expression for 
the spinless current (|in|) . which lead to corrections to the current which are of comparable magnitude 
as those which result from the wave focusing. Such terms are not explicitly present in the Dirac 
expression for the current (|12j) below. It is instructive to have a closer look at the structure of the 
electron current in a potential. A Gordon decomposition of the electron current 

f = e^-ff'^i' (12) 

can be performed by using the Dirac equation [i'j^(dp + ieAp) — m]^ = 0. The current can be split 
into a convective current and a spin current 

= _ (a^^)^]-^^'^A''| 

+ ^d,[^a^^^, (13) 

which are separately conserved and gauge invariant. For exact solutions of the Dirac equation, the 
two forms of the current are of course equivalent. Here, the convective current has exactly the same 
structure as the current of a scalar particle in a potential. Kopeliovich et al. calculate Coulomb 
corrections by describing the leptons as spinless particles. The discussion above clarifies when such 
an approximation is justified (see also sect. 3). 

In the following, we will consider electrons with energies of the order of several GeV. Then the 
impact of the focusing effect on the matrix element will not be larger than about one percent even 
for heavy nuclei like lead, where the electrostatic potential in the center of the nucleus is {Z = 82, 
A = 208) 

r.±'^ 25MeV, (14) 

2Ra ' ^ ' 

and Ra ~ 7.1 fm is the equivalent radius of a homogeneously charged sphere, which we will call 

"nuclear radius" for short in the following. E.g., for a scattering process where the initial and final 

momentum of the electron is of the order of 10 GeV/c, the focusing enters the cross section by a factor 

of the order {p' /pY = (10.025/10)'' ~ 1.01, according to (HT}. 
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3 Coherent Electroproduction of Vector Mesons 



3.1 Matrix Element 

Coherent electroproduction of vector mesons from virtual photons plays an important role in the 
understanding of the transition of soft diffr active models to quantum chromodynamics |121ll3j . Models 
based on the assumption of vector dominance have been applied successfully to describe the data Jl] . 
Our intention is not to give a better model to describe the vector meson production, but to explore 
the importance of Coulomb correction effects. Therefore we use a schematic model, which captures 
some essential features of the production process. But it is clear that for a more realistic analysis a 
better model should be used. 

Our model is the one proposed in 6 , which is inspired by the vector dominance model and which 
allows to derive a relatively simple form for the vector meson production amplitude on a nucleus with 
mass number A and nuclear charge Z. 

We denote the energy momentum vectors of the initial and final electron by (ei_2)i'i,2) and the 
momentum of the produced meson by pv- ey denotes the polarization vector of the meson. The 
production amplitude is then given by 

oo 

M{eA ^ e'VA) = j dxey ■ f{pi,P2,Pv,x), (15) 



where / = /i - /a and 



fl,2{Pl,P2,PV,x) 



_2 d_ 



X J eiP2-e2Pi + eiVx2(0 + e2Vxi(r) | 

X exp[i/?r + ixi{^ + iX2{r) — iijJi^2f]. (16) 

For details concerning the derivation of eq. (|lHjl we refer to P| and the appendix of this paper. We 
use a different sign convention for the eikonal phase than |S] . The gradient terms are artefacts of the 
spinless treatment of the electrons in ^ and therefore their physical significance in eq. dubious at 
best. However, from our discussion above follows that they can be neglected for momentum transfer 

= [pi — P2)^ — (ei — (-2) ^ -^A^' since they are of the order of uZ/Ra^ and \e1p2 — (-2Pi\ = V ^i^2Q^- 

In (|16() . we have 

ojI = {1- x)(ei - €2^ - x(l - x)p^ - 2x/B, 

ul = (1 - x)((ei - €2? - mlf - x{l - x)p^ - 2x/B (17) 

and K = pi — P2 — xpv- Note that in eq. (50) in [2], the squares for uJi^2 are missing and the terms 
containing B are wrong. We give therefore a detailed derivation of eqns. I17|) in the appendix, 
which is missing in jGj. Note also the different sign in front of uji^2 in the exponent. The misprints will 
be corrected in the electronic preprint version of the paper [Jj. 

B is the slope parameter of the differential cross section. For coherent electroproduction, it is 
related to the in68in chcirgc nuclca-r radius scjucircd. {f^j[)rms 

by (see e.g. [HI El) 

B = ^{rl)rms = lR\- (18) 
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At high energies and ei^2 S> \/Q^, the vectors pi, p2 and pi — P2 are nearly parallel. Therefore we 
choose the z-axis along pi such that the vector r = (6, z) is given by its z-component and the projection 
b to the normal plane. Using the relation 



oo 

dz 



r 

oo 



2E:o6J4-a;2, (19) 



where Kq is the modified Bessel function, and k,l,t the longitudinal and transverse components of k 
with respect to pi, we obtain {x = Xi + X2) 

00 27r 

h,2 = ~ , ^ fdbf dy,i^of6x/4^^)e--^™^(^)+^^(^) (20) 

LOl 2 C'^i^l 2 J J ^ ' ' 



-2 , ,2 



Eq. (|19|) is valid also for k| — a; < 0, whereas in |6j the sign in front of uj is problematic. From 

J dv.e^-^^™^(^) = 2^Jo(kt6), §;Ko{z) = -K^{z) (21) 





follows 



/i,2 = '"^:f "f^^ / dbbK,{b^4:^Unrb)e^-i^\ (22) 

4 - ^1,2 



Kopeliovich ei aZ. [H] performed calculations for a model potential 

VWo„o(r) = -— e-^'-(l - e-^^), (23) 
r 

with an infrared cutoff exp(— Ar), and the last factor corresponds to the monopole form of the nuclear 
form factor via 

^-2 ^ {r\)rras _ ^24) 

The eikonal phase can be calculated analytically for such a potential. One obtains 



X{r) = Xi{r) + X2{f) = - j dz Vmonoib, z) = 

— oo 

2aZ{i^o(A6) -i^o[(M + A)6]}. (25) 

We checked the result given in Fig. |1] in [Sj for small A. Our results show the same behavior, 
although we obtain slightly smaller Coulomb corrections. This might be due to the fact that the 
authors used formula (jlBf) only for their Sommerfeld-Maue calculations [7]. We can reproduce a 
nearly identical curve (Fig. if we reduce the slope parameter B according to formula ((TH|) by a 
factor of 2. 

For the nuclear radius Ra we used the formula 

Ra = (1.128 fm)yl^/3 + (2.24fm)^-^/^ (26) 
which is a good approximation for A > 20, whereas the mass number and charge are related by 

1.98 + 0.015^2/3- (27) 
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Figure 1: Coulomb correction for forward production for the model potential used in [Hj, calculated 
for ei = 100 GeV, y = {ei — e2)/ei = 0.6, and electron scattering angle 1° {Q^ ~ 1.2 GeV^), with the 
slope parameter replaced hy B ^ B /2 (see text). 

3.2 Improved Electrostatic Nuclear Potential 



It is obvious that the potential given by eq. (|23j) provides only an inaccurate description of the 
electrostatic nuclear potential near the nucleus. We looked therefore for potentials which allow an 
analytic calculation of the eikonal phase (x and xi,2 as well). This is not a trivial problem, since 
the class of meaningful analytic potentials which allow to express the eikonal phase by known special 
functions is rather restricted. 

A good choice is a potential energy for electrons of the form 

{aZ) Vmodelir) - - + ^2)3/2 " ^ (^2 + i?/2)4 ' (^8) 

which goes over into a Coulomb potential for r — > oo, and being close to the potential generated by 
the relatively homogeneous spherical charge distribution of a nucleus. The charge density 

p{r) = --d^irVmodelir)) (29) 
er 



corresponding to satisfies 



{p) = eZ, (30) 
{r^p) = ^R^eZ, (31) 
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i.e. we can identify with the equivalent radius of the homogeneously charged sphere which is given 
approximately by the formula given above. R' serves as an additional fit parameter. A good choice is 
R' = 0.5174i?. 

Because the eikonal phase turns out to be divergent for a Coulomb-like potential, we regularize 
the eikonal phase by subtracting a screening potential Vscr ~ (^^ + a?)~^/'^ from (|28|) . such that the 
potential falls off like for large r. The divergence can then be absorbed in a constant divergent 
phase ~ log(a) without physical significance, when the limit a — > oo is taken. It is instructive to 
calculate the eikonal phase for the simple screened potential 

V{r) = , r = + , (32) 

One obtains for a particle incident parallel to the z-axis with impact parameter h (r^ = 6^ + z^) 

1 1 



X? = aZ r dz 

J —oo 



aZ log , , (33) 

(z + Vr2 + a2)(62 + i?2) ^ ^ 



and therefore for the regularized eikonal phase 

■z + \/r2 + i?2. 



Xi 



^im (X? - aZlog(a)) = a^log( "y^^, ). (34) 



For a simple potential ~ (r2 + i?2)-i/2^ ^j^e rms radius does not exist, since the corresponding 
charge distribution does not fall off fast enough. 

For the potential Vmodei the regularized eikonal phase is given by 



f 7^-1 (t.^ 1 f R^ + b\ R^ 3 R''RiiRf + 4Rlb^ + 8b^) 

''^'^='<-^)-WT¥-Yo iRl + 62)7/2 (35) 



The z-dependent formulae for xi,2 are a bit lengthy, but can be derived in a straightforward manner. 

In Fig. [2 we compare the potentials generated by eq. the potential Vhcs of a homogeneously 
charged sphere with radius Ra and our model potential (|28j) for '^^^Pb with identical mean squared 
radii in all three cases. 



4 Results and Conclusions 

For the results presented in Figs. 3-5, we used the same kinematic conditions as for Fig. ^ Fig. 
El compares the eikonal correction, using the simple model potential Vmono and the potential Vmodei 
given by eq. ()28() . It turns out that the Coulomb corrections are overestimated by the model potential 
given by eq. (23). This is probably due to the fact that the potential is too deep in the central nuclear 
region. It is therefore mandatory to use an accurate description for the nuclear Coulomb potential in 
order to obtain reliable results for the Coulomb corrections. 

Fig. |1] shows the Coulomb corrections for each element with slope parameter B according to eq. 
(|18|) . but with three different rms charge radii for the electrostatic potential: the correct rms radius 
"v/I'^Drms, a too large rms radius 2^{r'jj^rms and a too small rms radius ^\/(r^)rms- The figure 
clearly indicates that the distortion of the electron waves is stronger for a small nucleus, whereas for a 
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Figure 2: Comparison of different model potentials used in the calculations. 



larger nucleus, the eikonal phase varies less over the length scale given by the nuclear radius (or slope 
parameter). The initial and final state wave function resembles then more a plane wave in the vicinity 
of the nucleus. The strong dependence of the Coulomb corrections on the size of the nuclear radius 
clearly indicates that the use of Sommerfeld-Maue wave functions would lead to incorrect results. 

The dependence of the results on the slope parameter is displayed in Fig. [SJfor '^^^Pb. There we 
varied the slope parameter B — > XB for A = 0.4 . . . 1.6, where B is the theoretical value given by eq. 
(|18jl. The results show a clear dependence of the Coulomb corrections on the model for the hadronic 
current. 

Finally, Fig. El shows the dependence of the Coulomb corrections for ^^O and '^^^Pb, where we 
have varied artificially the nuclear charge of the two elements, while the correct charge radius of the 
two nuclei and the corresponding slope parameter were held fixed. 

The calculations presented in this paper treat only the case where the scattering angle of the 
electron is small, such that the expression for the vector meson production amplitude can be reduced 
to a simple two-dimensional integral, which can be solved without involving large computational 
efforts. But it is also possible to perform numerical calculations of the three-dimensional integral 
representation of the amplitude on a modern workstation, such that arbitrary scattering angles and 
more general models of the hadronic current could be treated. This is the subject of a forthcoming 
paper. 

The oscillatory behavior of the corrected cross section shown in Fig. 3 in |S] can not be reproduced 
by our calculations. Fig. 4 shows clearly how important it is to use the correct charge distribution 
of the nucleus for the calculation of Coulomb corrections. Approximate wave functions for pointlike 
nuclei are therefore not adequate. 
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Figure 3: Comparison of Coulomb corrections for different potentials, given by eq. (|28)) . solid line, 
and eq. (PS)) . dashed line. Equal rms radii were used in both cases. 

5 Appendix A. Matrix Element for Electroproduction of Vector 
Mesons 

We start from the amplitude for coherent electroproduction of vector mesons by electrons derived in 
[S], eq. (27), working with the Born approximation first 

- ] ^^(Q^Py)- (36) 

where j is the spatial operator of the electron current, J the hadronic current operator, ly = ei — e2 
and Q'^ = — . From current conservation 

vk{q)=q-l{q) (37) 

we obtain 

M (x2 / nJ\.Q-,Pv)- (38) 

ei - £2 J q'^ - 

Using the hadronic model current (A = g — py here) 



m 



+ 1 + 5A2/2 



the amplitude becomes in Born approximation 

M oc 2(ll^^M)!V f ^3, f ^3 1 {2/BWy 

£1-62 7 7 {q^ -v^){q^ -v^ ^m\){2lB ^ 
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Figure 4: Coulomb corrections for the correct rms radius (solid line), doubled radius (dash-dotted), 
and radius divided by 2 (dashed). 

Going over to the eikonal approximation, we modify the electron current by the eikonal phases. 
Therefore, we have to evaluate the integral 

^ J J (g2-zy2)(^2_^2+^2^)(2/B + A2) ' ^ ' 

Due to the trivial identity 



my _ 1 



we can decompose integral I = Ii + I2 according to 



= I (fr I __Qiipi-P2-q)r+ixi(r)+iX2{r) ^ 



(42) 



(g2-z.2)(2/S + A2) 

In = f (fr ( g*(Pl-P2-^)r+ixi(r)+iX2(r) /^gx 

J J (^-zy2 + ^2^)(2/S + A2) ■ ^ ' 

With Feynman's trick 

L=j ^ (44) 



follows for h first 



g*(pi-P2-gOr+ixi(r)+«X2(r) 



h = j dx j d^r j d^q 



[(g2 _ i.2)(l _ ^) + + (g-- py)2)xl2 
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Figure 5: Coulomb corrections for different values XB of the slope parameter for 



pi(,pi-p2-q)r+ixi{r)+iX2if) 





Shifting the integration variable q ^ q + xpy leads to 



xpv) ^ — x^py ~ (1 ~ x)z^2 + 2x/B + xpy 



n2l2 ■ 



-p2 -xpv-q)T+ixi ir)+iX2 (r) 

h = I dx I d-'r I d-" 







[g2 - (1 - x)z^2 + _ a.)^2 _^ 2x/S] 



2 



ll" r n 1 d /"o e''^P^-P'2-^Pv-Q)r+iXi(r)+iX2(r) 



fdx f d'r-^ I d^q 



2 J J LOl duJi J (p — f^? 



Finally, the identity 



d q^ ^ = vr— 

q^ -ujf 



immediately gives the final result in agreement with eq. H16() 



1 



_ j dx ^ /" ^ ''' ^i(pi-P2-xpv)r-iuJir+ix(r) _ 

duJi J r 



For /2 we must simply replace ^ — niy. 
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Figure 6: Dependence of Coulomb corrections on the nuclear charge, but with approximately correct 
nuclear radii according to eq. and correct slope parameter for lead and oxygen. 
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